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OBJECT  OF  TASK 

To  improve  existing  knowledge  of  gamma  and  neutron-shielding  properties  of 
shelters,  and  where  necessary,  to  verify  experimental  ly  the  theoretical  information 
developed  in  this  field. 


ABSTRACT 

For  radiotion  shielding,  underground  or  buried  fallout  shelters  hove  an  Important 
advantage  over  other  types  of  shelters,  because  the  attenuation  of  the  radiotion  in 
such  o  shelter  is  primorily  a  function  of  the  thickness  of  the  material  in  the  roof  only. 

This  study  was  mode  to  investigote  the  dose  ottenuation  of  fallout  gamma 
radiation  by  various  thicknesses  of  concrete  roofs  of  buried  fallout  shelters  os  a 
function  of  time  after  a  nuclear  detonation.  A  spectrum  of  energies  is  used  for  the 
fallout  source  rather  than  a  single  average  energy  as  has  been  done  in  previous 
studies.  Dose  attenuation  factors  are  derived  and  presented  as  a  function  of  the 
above  parameters.  The  Office  of  Civil  and  Defense  Mobilizotion  recommends  a 
two-week  shelter-stay  time  in  the  event  of  a  nucleor  attock;  therefore,  also  presented 
is  an  average  dose  attenuation  factor  for  any  fourteen-day  stay  time  as  a  functian 
of  time  of  arrival  of  the  fallout  or  of  shelter-entry  time,  for  vorious  roof  thicknesses. 


INTRODUCTION 


In  this  nuclear  dge,  warfare  and  defense  problems  have  become  increasingly 
more  scientific.  A  nuclear  weapon  explosion  results  in  earth  and  bomb  debris 
which  is  contaminated  with  radioactive  fission  products.  This  radioactive  debris 
is  known  as  residual  radiation  or  fallout,  and  it  constitutes  o  serious  hazard  to 
unsheltered  personnel.  It  is  necessary  therefore  to  provide  shelter  from  the  harmful 
radiation. 

Underground  or  buried  shelters  have  a  definite  advantage  over  other  types  in 
resisting  the  effects  of  atomic  weapons,  especially  fallout.  The  amount  of  protec¬ 
tion  received  from  an  underground  shelter  is  a  function  of  the  moss  density  of  material 
in  the  shelter  roof  since  the  amount  of  radiation  coming  through  the  wolls  is  negli¬ 
gible  compared  to  that  coming  through  the  roof.  This  results  in  greater  net  protection 
than  afforded  by  a  surface  shelter  of  equal  cost. 

If  the  size  of  the  shelter  is  large,  that  Is,  not  the  small,  single-family  type, 
the  roof  can  be  approximated  by  an  infinite  concrete  slab  for  acceptably  precise 
mathematical  determination  of  its  shielding  capabilities. 

The  purpose  of  this  study  Is  to  investigate  the  dose  attenuation  foctor  for 
fallout  gamma  radiation  as  a  function  of  time  after  detonation  of  a  nuclear  weapon 
for  various  thicknesses  of  infinite  concrete  slab  shields.  In  addition,  using  the 
Office  of  Civil  and  Defense  Mobilization  criterion  of  a  two-week  stay  time  in  fall¬ 
out  shelters,  a  factor  has  been  calculated  which  will  determine  the  dose  received 
during  that  stay  time  as  a  function  of  the  H  +  1 -hour  dose  rote,  of  the  time  of  arrival 
of  the  fallout  after  detonation  or  of  shelter  entry  time,  and  of  the  slab  thickness  of 
the  shelter  roof.  An  average  dose  attenuation  foctor  for  any  fourteen-day  shelter- 
stay  time  as  a  function  of  the  above  Is  also  presented. 

Previous  work  has  been  done  by  Chilton  and  Saunders^  to  determine  the  dose 
attenuation  factor  as  n  function  of  roof  slab  thickness  by  using  an  average  energy 
of  1.0  mev  for  the  fallout  radiation.  This  study  will  to  a  great  extent  parallel  the 
above  work,  but  will  use  the  gamma  spectral  data  of  Nelms  and  Cooper^  for  various 
times  after  fission  to  specify  the  energy  of  the  fallout  rodiation. 


PROBLEM  CONSIDERATIONS 


The  geometrical  situations  investigated  are  shown  in  Figure  1.  In  Figure  1(d)/ 
it  is  assumed  the  fallout  is  evenly  distributed  on  top  of  a  smooth  infinite  plane 
surface.  Two  dose  points,  and  D2/  are  indicated.  D2  is  the  dose  received  at 
the  standard  3-foot  distance  above  a  uniformly  contaminated  infinite  plane.  D]  is 
the  dose  received  at  a  vertical  distance  h  beneath  the  contaminated  plane  with  a 
varying  thicl<ness  of  material,  t,  in  between  the  contaminated  plane  and  the  dose 
point  D].  In  this  study,  concrete  was  the  material  considered,  but  earth  could  be 
used  if  the  appropriate  mass  density  equivalent  is  used. 

It  has  been  found  that  in  the  computation  of  the  dose  at  D]  the  value  of  h 
does  not  greatly  affect  the  result  provided  h  is'less  than  a  mean  free  path  in  air  and 
t  is  about  0.  25  feet  of  concrete  or  greater.  The  area  of  interest  in  this  study  is  for 
roof  slabs  with* thicknesses  equal  to  or  greater  than  0.25  feet;  therefore,  for  conven¬ 
ience,  h  was  arbitrarily  chosen  to  be  3  feet.  This  choice  was  mode  so  that  when  the 
same  computation  is  made  for  t  =  0,  the  dose  calculated  at  D]  would  be  numerically 
valid  for  the  dose  received  at  D2> 

A  dose  attenuation  factor  for  the  smooth  plane  case  can  now  be  defined  os 
the  ratio  of  the  dose  received  In  the  open  at  a  distonce  of  3  feet  above  o  uniformly 
contaminated  plane  source  to  the  dose  received  inside  a  buried  shelter  where  the 
roof  approximates  a  concrete  slab  shield.  This  smooth  plane  dose  attenuation  factor, 
AF^,  is  equal  to  D]/D2* 

In  order  to  account  for  the  roughness  of  terroin  that  would  be  encountered  in 
practical  coses,  a  method^  suggested  by  the^U.  S.  Naval  Radiological  Defense 
Laboratory  •in  1955  is  used.  As  can  be  seen  in  Figure  1(b),  the  contamination  is 
assumed  to  be  uniformly  distributed  in  the  top  one-holf  inch  of  soil  under  the  infi¬ 
nite  plane.  The  distance  h  is  the  vertical  distance  from  the  incremental  segment 
of  contamination  to  the  dose-point  D3.  The  distance  M  is  the  distance  from  the  top 
of  the  soil  to  dose-point  D3  and  was  assumed  to  be  3  feet' in  the  calculation. 

Tf^e  dose  attenuation  factor  for  the  rough-plane  case,  AF2/  is  defined  as 
D]/D3.  (The  concrete  roof  of  the  structure  is  not  considered  to  be  rough,  even 
though  the  surrounding  plane  area  of  soil  is  so  considered.) 


CALCULATIONS 


The  calculaMons  used  in  this  study  for  the  doses  at  D|/  D2/  ond  D3  ore  slight 
modifications  of  the  calculations  presented  in  the  Chilton  and  Saunders  study. 
Details  of  the  derivations  are  presented  in  the  Appendix.  Only  the  three  integrated 
dose  equations  used  in  the  study  will  be  presented  here. 
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where:  Hi  =  Hq 


Ho  +  He- 


=  Photons/sec  for  each  energy  group  specified  by  Nelms  and  Cooper^ 

=  The  mean  energy  of  the  energy  group  at  which  ail  other  energy 
dependent  properties  are  evaluated,  mev 
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pQ  =  Linear  energy  absorption  coefficient  of  air,  cm'^ 
p  =  Density  of  air,  gm/cc 


p  =  Linear  total  absorption  coefficient  of  air 
(12.06  X  10”^  gm/cc),  cm”^ 

p  =  Linear  total  absorption  coefficient  of  earth  (1.442  gm/cc), 
6 

•  cm  ' 

Linear  total  absorption  coefficient  of  concrete 
(2.357  gm/cc),  cm“^ 


Exponential  integral  of  form  /  -  dt 

•  •  A 

•  A^,B^,A2/B2  =  Build-up  factor  coefficients  of  Berger  ond  Spencer. 

The  build-up  factor  is  defined  as  the  ratio  of  some  measurable  property  of  the 
photon  beam  (i.e.,  intensity,  number  of  photons,  energy  flux,  or  biological  dose), 
when  the  effects  of  all  quanta  are  included,  to  that  obtained  when  only  the  uncollided 
flux  is  considered.  Many  analytical  functions  have  been  derived  to  describe  the 
dose  build-up  factor  for  all  photon  energies.  In  this  study  a  biological  dose  build-up 
factor  of  the  following  form  (suggested  by  Berger  ond  Spencer^)  is  used: 

e 


t 

-Ei(-pX)  = 


-(B^-l)pr  -(B2-1)Hr 

B^  =  1  +  A]  B|  pr  e  +  A2  B2  pr  e  (4) 

where  A],A2»B],B2  are  dimensionless  coefficients  and  pr  is  the  number  of  mean 
free  paths  of  material.  The  coefficients  A^,B],A2,  and  B2  allow  the  empirical 
dose  build-up  factors  of  Berger^  to  fit  the  dose  build-up  factor  data  of  Goldstein 
and  Wilkins^  for  aluminum.  This  is  also  considered  reasonably  valid  for  concrete 
and  earth.  • 


RESULTS 

The  doses  from  the  above  equations  were  calculated  on  an  IBM-705  computer, 
using  the  spectral  data  of  Nelms  and  Cooper,  ^  the  air  linear  energy  absorption 
coefficients  provided  by  Berger,^  and  the  linear  total  absorption  coefficient  data  of 
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Gladys  W.  Grodsteln^  for  each  energy  group.  The  sum  of  the  doses  received  from 
each  of  the  energy  groups  represents  the  dose  received  through  a  particular  thickness 
of  slab  shield.  Fallout  spectra  for  1.12  hours,  5. 15  hours,  23.8  hours,  2. 13  days, 

4. 57  days,  9.82  days,  21 . 1  days,  45.3  days,  97.3  days,  and  208  days  after  fission 
were  investigated.  Data  for  these  spectro  are  tabulated  by  Nelms  and  Cooper.^ 


•  *  • 

Figure  2  is  a  presentation  of  the  dose  attenuation  factor  AFj  plotted  os  a 
function  of  concrete  shield  thickness  for  the  1.12-hour,  23.8-hour,  4. 57 -day, 

21 . 1-day,  and  208-day  spectra.  It  can  be  noted  that  the  attenuation  of  the  21 . 1-day 
spectrum  is  definitely  less  than  that  of  the  23.8-hour  spectrum.  On  the  other  hand, 
on  the  basis  of  data  provided  by  Miller,  based  on  a  spectrum  identical  to  Nelms  and 
Cooper^  and  plotted  in  Figure  3,  it  can  be  seen  that  the  average  energy  per  photon 
from  fission  product^  at  23.8  hours  is  greater  than  at  21 . 1  days.  If  the  mean  energy 
per  photon  is  a  good  estimation  of  the  penetration  power  of  the  radiation,  then^it  . 
would  be  expected  that  the  attenuation  of  the  23.8-hour  fallout  spectrum  would  be 
less  than  that  of  the  later  21 . 1  day  spectrum.  It  is  shown  in’Figure  2  that  the  oppo¬ 
site  is  true.  This  anomaly  is  discussed  in  the  next* section. 


Figure  4  shows  the  attenuation  factor  for  the  rough-plane  case,  AF2,*  plotted 
against  slab  thickness  for  the  same  spectra  as  AF^.  It  was  found  that  the  only  dif¬ 
ference  between  the  AF2  and  AF]  curves  is  thot  the  attenuation  for  the  rough-plane 
source  is  less  by  a  factor  of  about  1 .6  at  all  times.  • 


Figure  5  shows  the  attenuation  factor  VkF^  plotted  as  a  function  of  time  after 
fission  for  various  slab  thicknesses.  It  can  be  seen  that,  as  the  shield  thickens,  the 
attenuation  factor  varies  more  radically  with  time.  It  can  be  noted  that  the  second 
maximum  reached  at  about  ^0  hours  is  never  greater  than  the  initial  maximum  at 
1.12  hours,  for  any  of  the  thicknesses  specified.  Thus,  the  1.12  hour  fallout  spec¬ 
trum,  which  has  been  used  in  many  shielding  calculations  to  represent  the  fallout 
spectra  in  general,  is  still  a  conservative  basis  for  use. 

It  is  desirable  to  define  a  factor  F  as  a  function  of  the  time  of  arrival  of 
fallout  or  the  shelter  entry  time  (if  fallout  has*already  arrived)  for  various  slab 
thicknesses,  so  that  the  factor  F  when  multiplied  by  the  dose  rate  at  H  +  1  hour  will 
give  the  dose  received  by  sheltered  personnel  in  a  fourteen-day  stay-time  after  the 
burst.  This  factor  F  is  derived  by  integrating  the  attenuation  factor  multiplied  by  the 
t“1 .2  decay  scheme  over  various  fourteen-day  periods.  If  time  is  measured  in  hours, 
it  can  be  seen  that  the  fourteen-day  stay  dose  is  given  by: 


t,  +  336 

<Do*«)i4  days  =  ^o  ^ 

^  ti  •  • 


(5) 


% 


where:  =  H  1  hour  dose  rote  in  the  open  [at  reference  point  2*in  Figure  1(a)] 

•  •  * 

AF  =  Attenuation  factor  at  time  t  for  concrete  roof  thickness  ^  * 

T[Note:  AF(t,0)  =  l]  •  * 


t^  .=  Tim*e  of  arrival  of  fallout  or  start  of  14-day  stay-time  If  fallout  has* 
already  arrived,  hours  * 


Figure  6  gives  a  plot  of  thft  factor  F  as  a  function  of  t^ .  The  factor  does  not 
vary  in  a  simple  mathematical  way  as  a  function  of  shelter  entry  time  (t^).  This  con 
be  seen  Sy  the  variation  in  the  curves  plotted  for  the  2-  and  3-foot  slab  cases  after 
about  40  hours. 

An  average  attenuation  factor  can  be  plotted  for  any  fourteen-day  stay-time  as 
a  function  of  slab  thickness  and  time  of  arrival*  of  Rollout  or  shelter-entry  time.  This 
average  attenuation  factor  is  determined  by  dividing  the  F  factors  for  various  slab 
thicknesses  by  the  F  factor  for  a  0.0-foot  slab  thickness  at  a  particular  time.  The 
formula  is: 


ti  +336 


AF  (t,T)  dt 


AF 


ave 


(6) 


I 


t^  +  336* 


AF  (t,0)  dt 


Figure  7  Is  a  plot  of  AFgyg  as  a.function  of  t^  for  various  thicknesses  of 
concrete  roof  shields.  The  dose  that  will  be  received  in  the  open  during  any  14-day 
period  after  a  nuclear  explosion  can  be*calculated  by  the  following  formula: 


Dose  jgyj  (open)  o.2 


(t^ +336)0;  2  J 


(7) 


where  D  and  ti  are  as  defined  before.  The  dose  received  inside  the  shelter  for  a 
14-day  stay-time  would  be: 


•  • 


^**14  dajfs  (inside)  ^^1^14  days  (open)  ^^ave 


(8) 
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Figure  2.  Aftenuotion  factor  versus  concrete  thickness  for  smooth  isotropic  plane  source  for 
various  fission  spectra. 


Jic* 


Figure  7.  Average  attenuation  factor  versos  time  of  arrival  of  fallout  for  various 
concrete  thicknesses.  » 


DISCUSSION  AND  CONCLUSIONS 


In  Figure  2,  the  Chilton-Saunders  results  for  1-mev  photons  ore  plotted  os  a 
dotted  line  for  comparison  ^ith  the  present  results.  It  con  be  noted  that  the  slope 
of  the  Chilton-Sounders  curve  closely  follows  that  of  th*e  208-day  spectrum.  This 
indicates  that* the.  208-day  fallout  spectrum,  after  a  few  inches  of  concrete  has 
filtered  out  the  softer  gamma  rays,  has  penetration  properties  similar  to  monoenergetie 
gamma  radiation  wi4i  an  energy  of  1  mev.  The  1 . 12-hour  spectrum  is  obviously  more 
penetrating  than  a  l-me*v  photon  beam,  after  about  the  first-foot  of  penetration. 

• 

As  previously  noted,  the  attenuation  at  23.8  hours  is  greater  than  that  at  * 

21 . 1  days  even  though  the  mean  energy  per  photon  is  greater  at  23.8  hours.  This 
explained  as  follows:  In  the  spectral  dose  calculations  (Table  I),  it  con  be  seen  that 
the  contribution  of  the  spectrum  to  D]  changes  quite  radioally  fromi23.8  hoyrs  to 
21 . 1  days.  About  94.  2  percent  of  the  total  contribution  to  D]  for  the  21 . 1-day 
spectrum  came  from  radiatlon»in  the  1.47  -  2.95-mev  initial  energy  range,  while  for 
the  23.8-hour  spectrum,  only  79.3  percent  of  the  total  dose  comes  from  that  same 
initial  energy  rang&  and  a*greater  contribution  comes  from  photons  of  lower  Initial 
energies.  Thus,  even  though  the  average  energy  of  the  photons  may  be  greater^for 
the  23.8-hour  spectrum,  the  contribution  to  the  dose  shows  o  peak  at  a  higher  energy 
in  the  21 . 1-day  spectrum,  which  explains  the  lesser  attenuation  of  the  2K  I’-day 
spectrum.  After  21.1  days,  the  attenuation  factor  follows  the  expected  path  (Figure  5) 
but  does  not  peak  again  around  208  days  as  would  be  expected  from  the  Miller^  dato 
(Figure  3).  This  data  indicates  that  the  mean  energy  per  photon  is  obout  the  same  at 
208  days  as  it  is  for  23.8  hours.  The  explanation  is  similar  to  that  given  above. 

In  conclusion,  it  can  be  seen  that  at  no  initial  entry  time  after  1.12  hours  will 
a  person  in  a  shelter  for  fourteen  days  receive  a  greater  dose  thon  if  fallout  arrives 
at  or  before  1.12  hours  after  fission,  and  entry  time  is  at  1.12  hours.  This  is  obvious 
since  there  are  no  maxima  greater  than  at  1.12  hours  in  the  curves  of  FiQure  7  for 
any  of  the  thicknesses  indicated.  • 

If  however  an  average  attenuation  factor  for  a  particular  fourteen-day  stay-time 
is  needed  for  accurately  calculating  the  dose  to  be  received,  some  care  must  be  taken 
since  the  average  attenuation  factor  may  vary  as  much  as  a  factor  of  two  for  stays 
commencing  at  some  later  time  than  1.12  hours  after  fission.  , 
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Appendix 


DERIVATION  OF  INTEGRATED  DOSE  EQUATIONS  FOR  A 
CONCRETE  ROOF  SLAB  SHIELD  COVERED  WITH  FALLOUT 


•  • 


CALCULATION  OF  THE  DOSE  RECEIVED  AT  D^  [see  Figure  l(a)l  , 
The  dose  received  at  D]  from  on  infinite  contaminated  smooth  plane  is: 


= 


KE  P  . 
o  '^o 


Alt  p 


dA 


area  r 


(9) 


where:  E^,  p^  and  p  are  as  defined  for  Equation  1 

=  Dose  build-up  factor  *  * 

dA  =  An  incremental  area  on  the  surface  of  the  plane,  cm^ 

# 

<# 

♦  -1 

*  P]  =  Effective  linear  absorption  coefficient,^ cm 

.  • 

It  can  l?e  shown  by  similar  triTingles  that  dA  =  2ffrdr.  Thus,  Equation  7  becomes: 


GO 


D^  = 


2p 


L 


dr 

’  ®rT 


(10) 


The  effective  linear  adtorption  coefficient,  p^,  takes  into  consideration  the 
absorption  of  the  air  and  the  concrete  slob.  Therefore,  p  ^r  can  be  expressed  as: 


Mif  =  Ha(h-t)sec0  +  p^  t  sec0 


01) 


where  p^  and  p^  are  as  defined  for  Equation  1,  and  h,  t,  and  9  as  shown  in 
Figure  1(a).  * 


since:  r  =  h  seed 


then: 


HI  =  Ha  <He'Ha)~ 


*  '  *  *  4 

A  dose  build-up  factor  will  be  used  os  recommended  by  Berger  and  Spencer:^ 


Bj.  =  1  +  B^  H  ^r  e  +  A2  B2  M  e  • 

where:  A],  A2t  B]/  B2  ore  dimensionless  coefficients.  * 

e 

Substituting  Equation  13  into  Equation  10,  we  hove: 


6.  ^ 


O  “o 


f  a,-B,  H  re'*’'*'* 
r  ^ _  L  '  '  r 


f 


AjBjM,re  f  ‘ 


Looking  at  the  first  integral  of  Equation  14,. it  con  be  seen  that  it  is  in  the  * 
form  of  the  exponential  integral.  •  *  • 


r”.  -H,r  d(M,r) 

I  '  -p-77  =  -EiH 


In  computing  the  second  integral 


since:  • 


-BiM/  "^IHiT 

)  =  A^  e  (-B^  H^) 


=  -A,  B^  Ml  e  • 


•  • 


r  =  n 


00 


=  A^. 


(16) 


Similarly  for  the  third  integral:  , 


r  A  A  ■®2'''''dr  .  '“2^'' 

J  A2  B2  Mir  e  *7  =  ^2  e 

^r  =  h 


(17) 


Combining  Equations  15,  16,  and  17,*  we  have: 


K  E  Mrt 
o  "^O 


‘  D,  =  - 

•  '  2p 


■E|(-Mih)  +  e 


+  A2  * 


-B2Mih 


(1) 


CALCULATION  OF  THE  DOSE  RECEIVED  AT  D2  *  * 

•  * 

*  The  equation  for  the  dose  received  at  D2  is  derived  in  the  same  way  as  dose 
D^  except  that  only  the  absorption  by  the  air  need  be  considered.  Thus  pi  =  Pg. 

The  standard  3-foot  height  above  the  contaminated  plane  was  selected  (see  Figure  1), 
but  is  represented^in  the  equation  by  h.  • 


D2  = 


K  E  p 
o  o 

2p 


-BiMgh 

-Ej(-Mgh)  +  Ai  e  .  +  A2e 


V) 


CALCULATION  OF  THE  DOSE  RECEIVED  AT  D3 

*  The  dose  at  D3  is  due  to  the  contamination  mixed  uniformly  in  a  thin  layer  at 
the  surface  of  the  plane  as  shown  in  Figure  1(b).  This  method  is  used  to  approximate 
surface  roughness.  The  general  equation  for  the  dose  at  D3  is: 


D3  = 


o  ’^o 


/•  r  B.. 


^t  =  0.  wea 


(18) 


19 


where:  K,  p,  and  dA  are  as  defined  previously 

ff  dt,  and  T  are  as  shown  in  Figure  1(b) 
p  ^  =  An  effective  total  linear  absorption  coefficient 


From  Figure  1(b),  it  con  be  seen  that 


H  =  h  -  t 


•  r  =  h  sec0 


p  ^r  =  pg  t  SCC0  -  Pq  (h  -  t)  iec6, 


then: 


Ml  =  Mq  +  (Me  '  Mq) 


>9 

^here:  p^  ar\p  p^  are  os  defined  for  Equation  3  ^ 
fl  =  Angle  between  r  and  h  • 

Since:  dA  =  2ir  r  dr.  Equation  18  becomes 


•  p 

^  _  o* 
°3.-  2pT 


r  I  J  •.•■’•''4* 

'  '^t=0  "'r  =  h 


Using  the  build-up  factor.  Equation  \3,  end  substituting  into  Equation  20,  gives; 


l#'  c  -  *  ^ 

*  2pT  U,^o  '^r  =  h 


+  J  /  •  A,  B,  |i,re  ^dt 

t  =  0  r  =  h 

+  /  /  A2B2Pir  e  ^  ^  ^dt 

t  =  0  r  =  h  ^ 
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,  CotnputaHon  of  the  first  integral: 


f  f  =  f  -E,(-M,h)dt. 

■'  t  =  0  ^  r  =  h  ''t  =  0  • 


♦ 

Using  integration  by  parts  of  the  form:  fu  dv  =  uv  -  /v 

\J  ‘J 


I  * 

/_  -Ei(-Kl'')<f'  =  t 


0  t  =  0 


i  [E,(-H,h)]  d.  (23) 


p ,  h  =  H  +  p_  t 


J>lh)  =^p, 


d(Mih)  =  p^dt 


dt  = 


d(p,h) 


Also  when: 


t  =  0,  p^  h  =  p^H 
t  =  T,  p^  h  =  p^H  +  p^  T 


•*p^h 


Since: 


=  *f - d(-p  h) 


dr  #1  d 

~  Ej(-p,h)  =  - ! - (-P  h) 

•  dt  L  '  ’  J  d(-p  h)  dt  ’ 
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i 


To  simplify,  let  H  =  Pq  H  +  p^  T;  ond  since  T  =  (P3  H  -  p^  H)/pg, 
Equation  25  becomes: 

T  • 

J  -Ej(-p^h)dt  =  — •  I  -Ej(-P3H)P3H 

t=0  ^^e 


-p^H  -P3H 
“Ej{-PaH)  pg  H  -  e  +  e 


1  (26) 


Computation  of  the  second  integral: 


T  .  CO 


-B^p^r  r  -B^p^h 

,re  =  j  A,  e 


A,  B|  p  if  e  —  dt  =  /  A^  e  dt  (27) 

t=0"r=h  ^  "t=0 


But; 


p,  h  =  Pq  H  +  Pg  T 


*T 


Ai  e 


-B^p^h  r  -Bi(PgH  +  Pgt) 


-  I 


t  =  0 


A^  e 


dt 


=  A^  e 


-BiMa'’  r  -BiH,' 

/  e  dt 

^  t  =  0 


A^  e 


-BlMoH 


^-BiPe» 

h  He  • 


Ai  -BiPgH 
e 


®1  H, 


1  -e 


-BlHeT 


(28) 
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Since  T  =  P3H'’|jg  then  Equation  28  becomes: 


* 


-B1H3H 


) 


(29) 


Computation  of  the  third  integral: 

« 


Combining  Equations  26,  29,  and  30,  we  have  as  the  solution  fo  Equation  13: 


K  E  H,, 

o  ~o 
2p  He  T 


(-E.  (-H3H)  H3  H 


(-Ma”)  Mo  H 
♦ 


-m^h  -H3H- 
e  +  e 


# 


(3) 


•• 
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Mr,  Richard  Pork,  Notional  Academy  af  Scloncas,  2101  Constitution  Avonua,  Washington,  D.  C. 

Mr.  E.  E.  Sholowits,  Protoctivo  Construction,  GSA  Building,  19th  and  F  Street,  N.  W., 

Washington,  0.  C. 

Lt.  Col.  Russell  J.  Hutchinson,  052921,  Offico  of  Area  Enginaar,  Saudi  Arabia,  U.  S.  A.  Enginaar 
District,  Trons-East,  APO  616,  Now  York 

a 

Copt  A.  B.  Chilton,  CEC,  USN,  U.  S.  Naval  Civil  Enginaaring  Laboratory,  Port  Huanema,  Calif. 

Cdr.  C.  A.  Grubb,  CEC,  USN,  Public  Works  Canter,  Nevy  No.  128,  FPO,  Sen  Fronclsee 

Copt  W.  M.  McLallan  ,  CEC,  USN,  Public  Works  Office,  U.  S.  Naval  Base,  Charleston,  S.  C. 

Copt  L.  N.  Saunders,  Jr.,  CEC,  USN,  Bureau  of  Yards  and  Decks,  Cede  D>400,  Washington,  D.  C. 

Lcdr.  C.  Curlane,  CEC,  USN,  District  Public  Works  Office,  I4th  Naval  District,  Novy  No.  I2B, 

FPO,  Son  Francisco 

1  Cdr.  J.  F.  Clorfca,  CEC,  USN,  Area  Public  Works  OHica,  Chesapaake,  U.  S.  Naval  Waapens  Plant,  ^ 

Washington,  D.  C. 
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DUtrlkutlwi  LUt  (CMit'4) 


Mr.  0.  H.  AlbrtiKt,  Tk*  PwmaylvmU  St«t«  IMivarsIty,  CalUf*  af  EntlM«rii>t  Arciilfaatura» 
Unlaaralty  Park,  Paimayivania 

Cayt  J.  H.  Barkar,  Jr,  CCC,  USN,  U.  S.  Naval  Mlaaila  Canlar,  Paint  Mufw,  Calif. 

CJr.  W.  J.  Qirlatanaan,  CCC,  USN,  Bum  aw  af  VarJa  anJ  Daaka,  C»4*  D>400,  Waahlnftan,  0.  C.  * 

Cayt  J.  H.  LaflanJ,  Jr,  CEC.  USN,  Paarl  Harkar  Naval  SklyyarJ,  Navy  Na.  12S,  PPO,  San  Pranelaea 

Cayt  W.  A.  McManwa,  CEC,  USN,  U.  S.  Naval  Air  Statlan,  Narfalk,  Va. 

LcJr.  J.  D.  AnJrawa,  CEC,  USN,  Dafanaa  Atamlc  Swyyart  Agancy,  Waaklngtan,  D.  C. 

CJr.  0.  P.  Cunning,  CEC,  USN,  Olatrict  Pwfcilc  Warka  ONlaa,  4lk  Naval  Diatrlet,  Naval  Baaa, 
PkllaJalghla,  Pann. 

LcJr.  W.  J.  Prancy,  CEC,  USN,  Bwraaw  af  YarJa  anJ  Daeka,  CaJa  IVSOO,  Waaklnglan,  D.  C. 

Mr.  G.  L.  Arkuthnat,  Watarwaya  Eagarlmant  Statlan,  Paat  Offlea  Ban  631,  V'ckakwrg,  Miaa. 

LcJr.  H.  A.  Lacka,  CEC,  USN,  Caaig  SawJlay  D.  Bwtiar,  U.  S.  Marina  Carga,  Navy  Na.  161,  PPG, 

San  Pranelaea 

C6r.  H.  L.  Mwrgky,  CEC,  USNR,  PlalJ  CanwaanJ,  Dafanaa  Atamle  Sugpart  Aganey,  SanJia  Baaa, 
Al'kugwargwa ,  N.  M. 

Leir.  W.  H.  Sturman,  CEC,  USN,  San  Pranelaea  Naval  ShlgyarJ,  San  Pranelaea 

a 

LeJr.  W.  A.  Walla,  CEC,  USN,  U.  S.  Naval  Raaarva  Officara,  Training  Carga  Unit,  Un!varalty  af 
llllnala,  Urkane,  III.  ^ 

LeJr.  E.  M.  SawnJara,  CEC  USN,  Bwraaw  af  Yartia  anJ  Daeka,  Ca4a  D.MO,  Waakingtan,  0.  C» 

Lt.  C.  P.  Kriakankargar,  CEC,  USN,  U.S.  Naval  Amgklklawa  Canatrwctlan  Battallan  Twa,  PPO, 

New  Yark 

Lctir.  B.  S.  Marrill,  CEC,  USN,  U.  S.  Naval  Statlan.  Navy  Na.  I3S,  PPO,  New  Yark 

LaJr.  H.  W.  Stagkana,  CEQ  USN,  Diatrlet  PwUla  Warka  Offlaa,  13lk  Nava:  Diatrlet,  San  Brwna,  Calif. 

La6r.  R.  C.  Vanea,  CEC  USN,  U.  S.  Naval  Sakaal,  CEC  Offlaara,  Part  Hwanama,  Calif. 

LcJr.  W.  0.  Wllaan,  CEC  USN,.PwUlc  Wgrka  Cantar,  Navy  Na.  ?26.  PPO,  San  Pranalaca 
LeJr.  N.  W.  QaaMata,  CEC  USN,  Navy  Nwelaar  Pawar  Unit,  Pt.  Balvair,  Va. 

LeJr.  J.  C.  LaDawn.  CEC.  USN,  U.  S.  Naval  Sekaal.  CEC  OHlcara,  Part  Hwanaaie,  Calif. 

LtJG.  L.  K.  Donavan,  CEC,  USN,  Navy  Nwelaar  Pawar  Unit,  Pt.  Balvair,  Va. 

Offlea  af  ika  Cklaf  af  Enginaara,  Oagartawnt  af  ika  Army,  Waakingtan,  D.  C..  Attn;  ENOMC>ED 
CanwaanJlng  Officar,  U.  S.  Army  Ckamleal  Carga,  Raaaerck  A  Davalagmant  Caevaand,  Waakln^an,  D.  C. 
Offlear  In  Ckarga,  CECOS  (Attn:  ADCE  Cawraa),  Part  Hwanama,  Calif. 

Hdg.  U.  S.  Air  Parea,  Olractar  af  Reaaarck  and  Davalagmant,  DCS/D,  Waakingtan,  D.  C.  (Attn: 

Camkat  Camgananta  Divlalan) 

Cammandar,  Air  Parea  Sgaclal  Waagana  Cantar,  Kirtland  Air  Parea  Baaa,  Alkwgwargwa,  H.  M. 

Cammandar,  Plaid  Cammand,  Dafanaa  Atamlc  Swggart  Aganey,  Baa  5100,  Alkwgwargwa,  N.  M. 

Diraelar,  Qvll  Effacta  Taat  Orawg,  Atamlc  Energy  Cammlaalan,  (Attn:  Mr.  R.  L.  Carakla), 

Waakingtan,  D.  C. 

a 

Cklaf,  Dafanaa  Atamlc  Swggart  Agency,  Wa^lngtan,  D.  C. 


No*,  of 
copioo 


OUtrlbution  List  (Cont*4) 


OHico  of  Civil  and  Dofonoo  Mobil! totlon,  Attn:  Mr.  Bon  Toylor<  Bottio  Crook,  Mirh. 

U.  S.  Atomic  Energy  Commitsion,  Toebnicol  Inlormotion  Sorvico,  P.  0.  Bo>  67,  Oak  Rl4go,  Tonn. 

Commonding  Officor,  Ckomieol  Worforo  Lokorotorios,  Army  Otomicol  Contor.  M4. 

*  Lcdr.  C.  R.  Wkippio,  CEC,  USN,  NROTC  Unit,  Univoroity  of  lllinoio,  Urbono,  III. 

Hoodgwortort,  Field  Commend,  Oofonoo  Atomic  Support  Agency,  Sondio  Bote,  Albugworguo,  N.M. 

Office  of  tbo  Chief  of  Engineer*,  Ooportmont  of  iho  Army,  T-7,  Grevolly  Point,  Wothington,  0.  C 
Attn;  ENGNB 

Commending  Officer,  Engirtoor  Retoorck  &  Dovolepmont  Loborotorio*,  Fort  Bolveir,  Vo. 
Commanding  Officor  &  Oiroctor,  U.  S.  Novel  Rodielogicol  Oofon**  Loborotory,  Son  Francisco 
U.  S.  Army  Chomicol  Cantor,  Nuclear  Oofon**  Loborotory,  Edgovrood,  Md. 

Mr.  W.  R.  Perrot,  5113,  Sondio  Corporation.  Sandio  Bo**,  Albuguorquo,  N.  M. 

Dr.  N.  M.  Navnnork,  Civil  Enginooring  Holl,  Univorsity  of  Illinois,  Urbono,  III.  ^ 

Mr.  Konnoth  Kaplan,  Broadview  Rosoorch  Corporotioin,  1811  Trousdclo  Dr.,  Burlingomo,  Colif. 
Prof.  J.  Noils  Thompson,  Civil  Enginooring  Doportmont,  Univorsity  of  Tooos,  Austin,  Tos. 

Mr.  William  J.  Taylor,  Torminol  Ballistics  Laboratory,  Abordoon  Proving  Ground,  Abordoon 
Proving  Ground,  Md. 
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